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1 Introduction 

It is well-known that the behaviour for large t of solutions of the Cauchy problem 

Ut = ^u + uP, (1.1) 
m(x,0) = mo(x) (1.2) 

depends on the value of the exponent p of the nonlinearity. Let us first recall the critical 
value oi p = pp = l + ^/n called the Fujita exponent which borders the case of a finite- 
time blow-up for all positive solutions (for p ^ Pf) and the case of the existence of some 
global bounded positive solutions (if p > pp). It is also known that the Sobolev exponent 
Ps = is critical for the existence of positive steady states that is classical solutions 
ip G Co(M") of the elliptic equation Aip + ip^ = on M". Such solutions exist only if 
P ^ Ps (see e.g. [1], [8] ). Moreover, for p ^ ps there is a one parameter family of radial 
positive steady states ipk, k > 0, given by 

Mx) = kMk'^\^\), (1-3) 
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where tpi is the unique radial stationary solution with V'i(O) = 1, which is stricly de- 
creasing in |x| and satisfies ^/'i(lxl) -^0 as |a;| oo (see [17]). 
Another important exponent 



n - - 1 

^-^^ = 'a o I 1- ^""^ ^ 

n — 4 — l\/n — 1 

appeared for the first time in [lOj where the authors studied problems with the non 
linearities of the form j{u) = a(l + buy for some a,b > 0. It is also connected with 
a change in stability property of positive steady states defined in p.3p . Indeed, Gui 
et al. [9] proved that for p < pjl, all positive stationary solutions ipk are unstable in 
any reasonable sense, while for p ^ pjl they are "weakly asymptotically stable" in a 
weighted L°°— norm. Results on the asymptotic stability of zero solution to (ll.ll) - (ll.2p 
can be found in |16) and in the references given there. 

Let us recall that for p ^ pj^ the family of the positive equilibria ipk, k > 0, forms a 
simply ordered curve. Furthermore, this curve connects the trivial solution if — > and 
the singular steady state for — ^ cxd, which exists for p > Pst = in dimensions 

n ^ 3 and has the form Voo{x) = L\x\~^^'-^~'^^ with a suitably chosen constant L (see (12.1 p 
below). It is also known ([Ej) that ii ps ^ P < Pjl the graphs of the steady states 
Ipk, < k < oo, intersect the graph of Voo, whereas for p ^ Pjl have ipk < ^oo, 
{0 < k < oo). 

Our main goal in this note is to prove asymptotic stability of the singular stationary 
solution Voo in suitable weighted L^— spaces using estimates of a fundamental solution 
to a parabolic equation with singular coefficients [TTl [T3] . 



2 Results and comments 

It can be directly checked that for p > Pst 
singular stationary solution of the form 



= "/(n-2) and n ^ 3 equation (II. ip has the 



Voo{x) = L\x\ P~i = \^-—^(n-2-—-^^j \x\ p-i, (2.1) 

which plays the central role in this paper. 

In particular, problem p.ip -f |TT2l) with a nonnegative initial datum uq, which is 
bounded and below singular steady state Voo, has the global in time classical solution 
(see [IZl Th. 20.5 (i)] and [H, Th. 1.1]). Moreover, following Galaktionov & Vazquez 
[3, Th. 10.4 (ii)], we may generalize that result and prove that if ^ Uo{x) ^ Voo{x) and 
Uo{x) ^ Voo{x), then the limit function u{x,t) = limAr^oo w^(a;, t), where = u^{x,t) 
is the solution of the problem 

Ut = Au + u^, u{x,0) = mm{uo{x), N}, 
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solves (ll.ip and u{-,t) G L°°(M") for all t > 0. By those reasons, in the theorems 
below we always assume that u is the nonnegative solution to the initial value problem 
(ll.ip - (ll.2p with the initial datum Uq satisfying 

^ uo{x) ^ v^{x). (2.2) 

In order to show the asymptotic stability of the steady state v^c we linearize fll.ip 
around v^o- Denoting by m = u{x,t) the nonnegative solution to (ll.ip - fll.2p and intro- 
ducing w = Voo — u, we obtain 

Wt = Aw + j^w - [{v^ - wY -vP^+ pyP^'w] , (2.3) 

where 

A = \in,p) = (n-2- -^). (2.4) 

p — l\ p — 1/ 

Next, we use estimates of the fundamental solution of the linear heat equation with 
singular potential 

ut = Au + j^u, X G M", t > (2.5) 

obtained recently by Liskevich & Sobol [11], Milman & Semenov [I3j (see also Moschini 
& Tesei [l2]). As the consequence of the Hardy inequality, it is crucial in that reasoning 
to assume that A ^ in equation f l2.5p . Coming back to the perturbed equation 

(12. 3p and using the explicit form of A(n,p) in (12. 4p . we obtain by direct calculation (see 
Remark [4?T] for more details) that the inequality \{n,p) ^ IZiziL jg valid if 



n — 2\/n — 1 

P>Pjl = ^ , ^ lor n ^ 11. (2.6) 

— 4 — — 1 

By this reason, we limit ourselves to the exponent p of the nonlinearity in (11. ip satisfying 
(12. 6p . The exponents mentioned above are ordered as follows: pp < Pst < Ps < Pjl- 

We introduce the parameter a which plays a crucial role in our reasoning by the 
formula 

. , n-2 I in - 2)2 2p ( 2 A , , 

(T = (T(n,p) = —in-2 . 2.7 

2 Y 4 p-l\ p-lj ^ > 

It is worth pointing out that a{n,p) > "^/(p-i) if p > Pst and n > 2. Moreover, the 
number cr{n,p) has the property 2a{n,p) < n. Let us also notice that cr{n,p) appears 
in a hidden way in the papers of Polacik, Yanagida, Fila, Winkler (see e.g^.f!!], [6j), 
because it is the sum of the constant y(p-i) and Ai, where Ai is one of the root of the 
quadratic polynomial — {n — 2 — 2L)z + 2{n — 2 — L), given explicitly by the formula 

\i = ]^(^-2-2L- V(n - 2 - 2L)2 - 8(n - 2 - L)^ , 

where L is defined in (12. ip . 

Now we are in a position to formulate our first result on the convergence of the 
solutions towards the singular steady state. 
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Theorem 2.1. Assume (12 ■2p . (12 .Gp . (I2.7p . Suppose, moreover, that there exist constants 
b > and i E (a, n — a) such that 



Voo{x) — b\x\ ^ Uo{x) 



for all \x\ ^ 1. Then 



sup \x 

\x\^Vt 



{voo{x) - u{x,t)) ^ Ct 



2 



(2.8) 



and 



sup (yoo{x) — u{x,t)^ ^ Ct 2. 



(2.9) 



for a constant C > and allt^l. 

Polacik & Yanagida [HI Th. 6.1] showed that under the assumptions of Theorem 
12.11 the pointwise convergence holds true, namely, limt^oo ^(a;, t) = Voo{x) for every 
a; e M" \ {0}. More recently, Fila & Winkler [5] proved the uniform convergence of 
solutions u = u{x,t) toward a singular steady state on M" \ -Bjy(O), where -Bj/(0) is the 
ball in R" with the center at the origin and radius u. Theorem 12. Ib ompletes those results 
by providing optimal weighted decay estimates in the whole M". 

Remark 2.1. Note that our calculations in the proof of Theorem 12.11 are valid for any 
i G (2/(p-i), n — a), but for i G (2/{p-i), cr] the right-hand side of inequality p. 81) does not 
decay in time. □ 

We can improve Theorem 12.11 for £ = a as follows. 

Theorem 2.2. Assume that (12. 2p . (12. 6p and (12.71) are satisfied. Suppose that there 
exists a constant b > such that 



Voo{x) - b\x\ " ^ Ui){x). 



Let, moreover. 



lim \x\'^ ivooix) — uq{x)) = (). 

\x\^oo 



Then 



lim sup \x 



''{v^{x)-u{x,t)) =0 



and 



lim t2 sup {v^{x) — u{x,t)) = 



\x\^Vt 
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Corollary 2.1. Under the assumptions of Theorem \2.1\ and Theorem \2.^ respectively 
if, moreover, b is sufficiently small, we obtain 

||M(-,t)||oo > Ct-(v~r)-2 £^(^(j^n-a) (2.10) 

for a constant C > and allt^l and 

lim ||M(-,t)||oo = +00 tf i = (7. (2.11) 

t—*oo 

Remark 2.2. Estimates from below of ||m(-, t)||oo, similar to that stated in (12. 101) . were 
obtained by Fila et al. in j2], Theorem 1.1.], [HI Theorem 1.1.] and improved in [H] 
Theorem 1.1.] using matched asymptotics expansions. In Corollary 12.11 we emphasize 
that this inequality is an immediate consequence of Theorem 12. 1[ □ 

Remark 2.3. For p > pjl estimates (12. 8p and (12.91) seem to be optimal, because they 
imply the optimal lower bound (I2.10p . see [3]. On the other hand, for p = Pjl the 
authors of [4] obtained the logarithmic factor on the right-hand side of ( 12.101) . which we 
are not able to see by our method. □ 

Our next goal is to prove the asymptotic stability of foo in the Lebesgue space L^(R"). 

Theorem 2.3. Assume that (Q, (I2l6|) and are valid. 

i) Suppose that - Uq e ^^(IR") and \-\'"{voo- Uq) G /.^(IR"). Then 

\M)-u{-,t)\\2 ^ Ct-'^\\voo-uo\\i + Ct~'^\\\ ■ r"(t;oo -^^o)||i. (2.12) 
a) Suppose that v^c — Uq E L^(]R"'). Then 

lim \\vooi-) - u{-,t)\\2 = 0. 

t—^O0 

Note, that Voo G L^^^(M") for every p > pp. Here, this property of the singular 
solution foo is satisfied, because pji > pp. 

Using the fact that the steady states ipk defined in (11.31) are below the singular 
stationary solution Voo for p > pjl, we may rephrase Theorem 12.31 as follows. 

Corollary 2.2. Assume that (12.21) . ( 12.61) and (12.71) are valid. Let ipk be the stationary 
solutions (|1.3p for some k > 0. 

t) Suppose that - Uq e L^{W) and \ ■ {'"{ipk - Uo) G ^^(R"). Then 

\\M-)-<,t)\\2 < crt||^fc-Mol|i + cr'^|| | • |-'^(^fc - no)||^. (2.13) 

a) Suppose that ipk — uq E L^(M"). Then 

lim ||^,(.)-«(.,t)||2 = 0. (2.14) 

t—>oo 

Remark 2.4. Observe that Theorem 12.31 and Corollary 12.21 complete the results by 
Polacik & Yanagida, who proved in [HI Proposition 3.5] the stability estimate 

\\i>ki-) - u{-,t)\\2 ^ \\ll>k - Molh- 

□ 
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3 Linear equation with a singular potential 

In this section we recall the estimate from above of the fundamental solution of the 
equation Ut = Au + X\x\~^u obtained by Liskevich & Sobol in [11] and by Milman & 
Semenov [I3]. Following those arguments, we define the weights (fa{x,t) e C(R" \ {0}) 
as 

, , f i^Y if |x| ^ v^, 

II II \x\ ^ y't. 

Theorem 3.1. [HI [I3] Let Hu = Au + X\x\-^u. Assume that ^ A ^ ("-2)74. The 
semigroup e^^^ of the linear operators generated by H can he written as the integral 
operator with a kernel e~^^{x,y), namely 

e-^"uoix)= e-^"{x,y)uo{y) dy. 

Moreover, there exist positive constants C > and c> 1, such that for allt>0 and all 
x,y eW\ {0} 

^ e-*^(x, y) ^ Cifi^ix, t) ^„{y, t) G{x - y, ct), (3.2) 



where a = — y ^^^^^—^ A, the functions ip^ are defined in (13.11) (see also Remark 

1/2 



\3.1\ below) and G{x,t) = (47rt) exp(-l^iy47rt) is the heat kernel. 

Remark 3.1. In fact, Milman & Semenov in [13] used the more regular weight functions 
e C2(M"\{0}), namely 

if |x| ^ 2^/t 



and I ^ $cr(a^,t) ^ 1 for ^ \x\ ^ 2\^. It can be checked directly that there exist 
positive constants c and C for which the inequalities 

hold true, where ipo- are defined by (13. ip . By this reason we are allowed to use the 
weights (fa instead of □ 

The following theorem is the consequence of the estimates stated in (13.21) . 

Theorem 3.2. Let the assumptions of Theorem \3.1\ he valid. Assume that p > 1 + 
Suppose that there exist b > and i G (^^, n — a) such that a nonnegative function wq 
satisfies 

2 

■^0(2;) ^ b\x\ p-i for \x\ ^ 1, 
""^0(2;) ^ for \x\ ^ 1. 

Then 

sup ip-\x,t)\e-^"wo{x)\ ^ Cr^ (3.3) 
for a constant C > and allt^l. 
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Proof. First, for every fixed x G M", we apply the estimate of the kernel e *^ from 
Theorem 13.11 in the following way 



Next, we split the integral on the right-hand side into three parts /i(x,t), l2{x,t) and 
l3{x, t) according to the definition of the weights ipa- and the assumptions on the function 
Wq. Let us begin with Ji(x,t): 



Ii{x,t)=C / G{x -y,ct)ipa{y,t)wo{y) dy 



71 — a 
2 



< Cbt^ / G{x-y,ct)\y\'''~^ dy ^ Cbt 

n — 2 

because G{x — y,ct) is bounded by Gt~2 and the function \y\ " p-^ is integrable for 

||/| ^ 1 if P > 1 + "^/{n-a). 

We use the same argument to deal with 



hi^x, t) = G G{x - y, ct)ip„{y, t)wo{y) dy 



cr — n 



^Cbt^ / G{x - y,ct)\y\~''-' dy ^ Gbt— / \y\ dy 

^ Gbr-2 +Gbt-^. 

Finally, we estimate 

hix, t) =G G{x - y, ct)ip„{y, t)wo{y) dy 



2 



'\y\^Vt 

^Gb [ G{x-y,ct)\y\-^ dy ^Gbt- 
J\y\^Vt 

using the inequality 1 ^ fo'^ \y\ ^ identity G{x — y, ct) dy = 1 for 

t > 0, X e M". Since ^ G {^'/{p-i), n- a), we complete the proof of (l331) . □ 

Theorem 3.3. Assume that \ ■ {"wq e L°^(]R") and 

lim \x\'^wq{x) = 0. 

|2;|— >oo 

Then 

lim sup (p~\x,t)\e~^^wo{x)\ = 0. (3.4) 
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Proof. For every fixed x G we use the estimate from Theorem 13.11 as follows 

(p~^{x,t)\e~*"wo{x)\ G{x - y,ct)ip^{y,t)wo{y) dy. 



We decompose the integral on the right-hand side according to the definition of ipa- and 
we estimate each term separable. Substituting y = we obtain 

h{x,t)^C [ G{x-y,ct)(^) wo{y)dy 



Cr^ I g(^- z,c\\z\-^''\Vtz\"wQ{Vtz) dz. 



'\z\ 

Hence, 

1 2 sup Ji (a;, t) as t ^ oo 

by the Lebesgue dominated convergence theorem, because G[-^ — z,c) is bounded and 
the function \z\~'^'' is integrable for \z\ ^ 1. By the assumption imposed on Wq, given 
£ > we may choose t so large that 

sup lyl^woiy) < e. 

Now, using the inequality 1 ^ {'^Y ^'^^ 1^1 ^ obtain 

l2{x,t) = G{x - y,ct)wo{y) dy ^r'^ G{x - y,ct)\y\''wo{y) dy 



^ et ^ G{x — y, ct) dy. 

J\y\^Vt 

Since Jjg,„ G{x — y, ct) dy = 1 for alH > 0, x G M" and since 5 > is arbitrary, we get 

t2 sup l2{x,t) — i> as t ^ oo. 

□ 

Let us defined the weighted L^-norm as follows 



and 



\f{x)(p^^{x,t)\'^(p1{x,t) dx ) for every 1 ^ g < cxd. 



.,v<.W = sup (^^^(a;,t)|/(x)| for q = oo. 
Note, that in particular for q = 2, the norm || ■ ||2,i/3^(t) agrees with the usual L^-norm on 



Proposition 3.1. Suppose that 1 ^ g ^ oo. Then the following inequality holds true 

||e"*'^wo||<?,vp<.(i) ^ Ct'^^~'^''\\wo\\r,<p^{t) (3.5) 
for every 1 ^ r ^ g ^ oo and all t > 0. 

Proof. The proof of estimate (I3.5p can be directly deduced from the reasoning by Milman 
& Semenov. Indeed, in [I3l page 381], we can find the inequality 

Hence, substituting (p^rf = Wq and using the definitions of the norm || ■ ||g,,^^(t), we obtain 
(13 -Sp with q = 2 and r = 1. This inequality together with 

stated in [131 P^^ge 391], imply (13 .51) for g = 1 and every 1 ^ r ^ 2 by Riesz-Thorin 
interpolation theorem. Moreover, the operator e~*^ is self-adjoint, so by duality the 
inequality 

1 1 tFf 1 1 — 1 1 II 

||e Wo\\ao,^^{t) ^ Ct ^\\Wo\\2,^„{t) 

holds true. The semigroup property e~*^ = e^'/^^e^*/^^ leads to (13. 5p with q = oo and 
r = 1. Applying duality and Riesz-Thorin interpolation theorem once more, we complete 
the proof of (13.51) . Let us emphasize at the end of this reasoning, that the inequalities 
( 13. sp are used by Milman & Semenov in [1^ to derive the kernel estimate (13. 2p . □ 

4 Linearization around a singular steady state 

Let M be a solution of (11.10 with initial datum satisfying (12.21) . We substitute 

w{x, t) = Voo{x) — u{x, t) 

to get 

wt = Aw + —w - [{voo - wY -v^+ pyP^^w] , (4.1) 

where A = X{n,p) = -^{n — 2 — ^^)- Let us note that the last term on the right-hand 
side of equation (14.11) is non positive, namely 

which is the direct consequence of the convexity of the function /(s) = s^. Indeed, since 
the graph of the function / lies above all of its tangents, we have f{s—h)—f{s) ^ —f'{s)h 
for all s and /i in M. 

The proofs of our results are based on the following elementary observation. If w is 
a nonnegative solution of equation (14.11) with the initial condition Wo{x) ^ 0, then 

^ w{x,t) < e'^^woix) 
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with Hw = Aw + X{n,p)\x\ "^w. Consequently, using the condition ^ Uq{x) ^ Voo{,x) 
and the just-mentioned comparison principle we can write 

^ v^{x)-u{x,t) ^ e-'^ {v^{x) - Uo{x)) (4.2) 

or , equivalently, 



Vr^\x] - e 



{voo{x) -Uo{x)) ^ u{x,t) ^ Voo{x). (4.3) 



Remark 4.1. If n ^ 11 and either p ^ pjr or < p < "t^^TV^ . then the linearized 
problem 

w{x, 0) = wo(a;), 

with A = X{n,p) = -^ziin — 2 — ^^) has the unique solution. Indeed, in the view of 
Theorem 13.11 it is sufficient to show that 

A(n,p) = — ^ n-2 -U^ ' 



p — 1\ p — 1 J 4 

Substituting y = V(p-i), after elementary calculations, we arrive at the inequality 

IQy^ + (32 - 8n)y + - 12n + 20 ^ 

which has the solution ye ( - oo, '^-^-^v^ y ^ Moreover, if n ^ 11, 

then "-4-2v^ > and if n G (2, 10), then »-4-2v^ < g and "-4+2v^ > q. These 
observations give us that p ^ Pjl or < p ^ "l^^w^^ - D 

5 Asymptotic stability of steady states 

Proof of Theorem \2.1i It suffices to use inequality (14.21) and to estimate its right-hand 
side by Theorem 13.21 □ 



Proof of Theorem \2.^ As in the proof of Theorem l2. 11 it is sufficient to use (14.21) together 
with Theorem 13.31 substituting wo{x) = Vooix) — uo{x). □ 
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Proof of Corollary \2.1\ Since we have inequality (I4.3p . it suffices to prove that 

sup [voc{x) - e~*-^wo(x)] ^ C(6)t-(f-i)-2 

for t^o = ?^oo — ^0- Hence, inequality (|3.3I) from Theorem 13.21 enables us to write 



for all a; G M" \ {0} and t > 0. Next, using the explicit form of the weights ipa-, we define 
the function 



F{\x\,t) = Vooi\x\) - Cbip„{x,t)t-2 



L\x\ p-i — for |x| ^ vt, 

— Cbt~^ for |x| ^ y/i. 



An easy computation shows that the function F has its maximum at 

|x| = C{b)t 2 <T(p-i)-2 

and it is equal to 

maxF(|x|,t) = C{b)t<^=^ 

for some constant C{b) ^0. Hence, we get (12. 101) . 

To obtain (I2.11I) . we use the result from Theorem I3.3[ It follows from (13 ■4p that for 
every e > there exists T > such that 

\e~^^wo{x) \ < eipcr{x,t)t~^ 

for all X eW'\ {0} and t > T. Hence, by (HJD, we have 

Voo{x) - e"*^(foo(x) - Uo{x)) ^ v^{x) - Ceip„{x,t)t~'^ . 

Now, once more using the explicit form of the weights (fa, we consider the function 

G{\x\,t) = Voo{\x\) - Cbipa{x,t)t~'^ 

An elementary computations give us that the function G attains its maximum at 

p-i 

\x\ = ce <^(p-i)-2 

and 

maxG(\x\,t) = Ce'-'^p-^^-^ 

for some constant C ^ 0. Since a > y(p-i), we see that the maximum of the function G 
diverges to infinity if e tends to zero. This completes the proof of (12.111) . □ 





\x\ 


1^1 


x\ 
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Proof of Theorem 1 2. 3\( i) . According to (I4.2p it is enough to estimate the L^— norm 
of the expression e~^^WQ for every Wq satisfying two conditions: Wq G L^(]R") and 
I ■ |~°^ifo G L^(W^). Applying (13 .SI) , with g = 2, r = 1 and using the definition of the 
functions ipa{x,t), we may write 



n — 2a 



||e wo\\2^Ct ''\\wo\\i,^^{t) = Ct 4 / wo{x)\x\ dx 

J\x\<:Vt 

n I n — 2a . i -r 1 1 1 1 1 1 

+ Ct"4 / wo{x)dx^Ct ^||wo| ■ r 111 + C'^"^||wo||i. 
This estabhshes formula fl212D . □ 



Proof of Theorem \2.3\( ii). Again, by (I4.2p . we only need to show that 

lim ||e"*'^M;o||2 = 



t^oo 



for each Wq G L^(M"). Hence, for every e > we choose ip ^ C^(M") such that 
11^0 ~ ^Ih < Using first the triangle inequality and next (13 .Sp . with q = 2 and r = 2, 
we obtain 

We-'^'woh < ||e-*^K - V^)||2 + ||e-*^V^||2 



^ C£+ lle-*^^ 



2- 



Since the second term on the right-hand side convergence to zero as t — oo by the first 
part of Theorem 12.31 we get 

lim sup ||e^*''^Wo||2 ^ Ce. 
This completes the proof of Theorem 12.31 (ii), because e > can be arbitrary small. □ 



Proof of Corollary \2.2l We linearize equation p.ip around the positive steady state ipk 
substituting v = ipk — u to get 

Vt = Av +pipl~^v - [{ipk - vY - ipl + pipl'-^v). (5.1) 

Once more, using the convexity of the function /(s) = s^, let us notice that the expression 
{ipk —vY — ipk~'^v is nonnegative. Furthermore, ipk < foo as long as p > pjl and n ^ 11. 
Hence applying first, the comparison principle to the approximate problem 

v{x, 0) = Vq{x) 

with the nonnegative initial datum and next, passing to the limit when N tends to 
infinity, we get 

^ ilJk{x) - u{x, t) ^ e~*^ {iJk{x) - uq{x)) . 

Now, (12.131) and (I2.14p are the straightforward consequences of the reasoning used in 
the proof of Theorem 12.31 □ 
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